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Experimental study of work exchange with a granular gas: 
the viewpoint of the Fluctuation Theorem. 
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Abstract - This article reports on an experimental study of the fluctuations of energy flux 
between a granular gas and a small driven harmonic oscillator. The DC-motor driving this system 
is used simultaneously as actuator and probe. The statistics of work fluctuations at controlled 
forcing, between the motor and the gas are examined from the viewpoint of the Fluctuation 
Theorem. A characteristic energy Ec of the granular gas, is obtained from this relation between 
the probabilities of an event and its reversal. 
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Introduction. — Often, the principle of probing is to 
measure the response of a system to an excitation imposed 
If om outside. This excitation must be small enough not to 
perturb the state of the system under scrutiny. Together 
vith the Fluctuation Theorem (FT), this principle is used 
in the present study to probe the disordered motion of a 
granular gas. 

In addition to the power needed to keep the granular gas 
ip a Non Equilibrium Steady State (NESS), (i.e. off- 
setting the mean dissipation), a small probing power is 
ijnposed. The response to this excitation is characteristic 
of the granular gas disordered state. More precisely here, 
i}; is the product of the forcing perturbation and the ve- 
locity response, i.e. the power perturbation, that carries 
ipformation on the system: its statistical properties will 
be considered in order to characterise the disordered gran- 
ijlar gas NESS. 

The measurement (excitation-response), is performed by 
means of a driven harmonic oscillator coupled to the gran- 
ular gas. The Fluctuation Theorem (FT) is used to com- 
pare the work given and received from the granular gas by 
the harmonic oscillator. It is the Steady State Fluctuation 
Theorem (SSFT) that is considered all along this article, a 
relation holding in the limit of asymptotically large times. 
No further mention of this distinction will be given later. 
The power provided (work flux) can be related to the rate 
of entropy production a. The FT states a relation be- 
tween the probability of events of positive a, and that of 



equal but negative entropy production rate. It states that 
the ratio of these probabilities simply increases as the ex- 
ponential of cr [T]. It therefore quantifies the failure of 
detailed balance, for a large number of degrees of freedom 
chaotic dynamical system. 

This experimental study applies an original and very sim- 
ple probing principle, to show that FT seems to holds in a 
stationary granular gas, and measure a characteristic en- 
ergy Ec. It is a contribution to the very few experimental 
studies of the different forms of the FT. (For a nice review, 
see H].) 

Rather few experiments invoking the FT have been actu- 
ally done in granular gases particularly. Most studies are 
numerical [3] , or theoretical [1] . Strictly speaking, the only 
experimental study of a granular gas with the FT is that 
of Feitosa and Menon [5] . In a 2D vibrated granular mate- 
rial, they measured by video tracking the momentum flux 
of particles in and out of a sub- volume, and used the FT 
to define an effective temperature of the medium. The 
interpretation of these results is difficult, and still ques- 
tioned [4]. 

In a distinct context, D'Anna et. al. performed rheol- 
ogy in dense granular fluids using a driven torsion pen- 
dulum. Thanks to an out-of-equilibrium extension of the 
Fluctuation-Dissipation Theorem, they defined an effec- 
tive temperature [6]. 

The present article describes an original study of a dilute 
3D granular gas, also using a driven torsion pendulum. In 
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a simpler but reliable manner, it show how to measure an 
energy Ec, characteristic of the particles disordered mo- 
tion, thanks to the FT. 

The equation of motion of the harmonic oscillator writes 
as the following: 



MO + -fO + ke = Tit) + riit), 



(1) 



where 9 is the angle of torsion of the pendulum, and dots 
stand for time derivatives. M, 7, k are respectively mo- 
ment of inertia, viscous friction coefficient, and spring con- 
stant. A sine torque T(t) is imposed from outside. The 
last term r]{t) represents the coupling with the NESS gran- 
ular gas heat bath. In the framework of Langevin equation 
description, it represents the momentum transfer rate at 
each shock from the beads. 

When the FT is expressed in terms of the work varia- 
tion rate w{t) ~ Td, the mechanical power transmitted 
during a time-lag t to the gas is expressed as: Writ) = 
- J w{t')dt' . It relates the probabilities of giving a power 
Wr to the gas and the probability of receiving the same 
amount from it. It states that this ratio increases expo- 
nentially with the coarse-grained work riiir ■ 
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for asymptotically large times r. 11 is the probability, and 
the coefficient Ec is a characteristic energy (possibly linked 
to the mean kinetic energy of the gas) . 
The principle of the measurement and the experimental 
set-up are described in the next section. The results ob- 
tained are detailed in the following section. The last sec- 
tion is devoted to a brief discussion of these results. 

Measurement's principle. — Fig. [1] sketches the ex- 
perimental set-up. The granular gas is composed of a few 
hundreds of 3 mm diameter stainless steel beads, contained 
in a vibrating vessel. 

The vessel is aluminum made, 5 cm diameter, 6 cm deep, 
its inside bottom is slightly cone-shaped to favour hori- 
zontal momentum transfer (angle=10°). (See fig. [H) The 
electromechanical shaker is driven by a sine generator, via 
a Kepco current amplifier. An accelerometer fixed on the 
vibrating vessel measures the vertical acceleration: be- 
tween 41 ms~^ and 60ms~^, at a frequency /oxc. = 40 Hz. 
A small DC-motor is simultaneously used as actuator and 
sensor. It is a regular permanent magnet, brushed DC- 
motor, of relatively small size (25 mm diameter). A plastic 
blade of approximately 20 x 20 mm is fixed on the axis of 
the motor. A torsion spring is used to produce an elas- 
tic force on the motor axis. This system motor + blade 
-f- spring forms a harmonic oscillator. Its resonance fre- 
quency is a few hundreds of Hz, in any case higher than 
any frequency of the signal of interest. The motor is fixed 
on a cover closing the vessel, which prevents the beads 
from hopping off the vessel. 
The principle of the measurement is the following. A DC- 
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Fig. 1: The mechanical system is composed of a vibrating vessel 
with its driving, and the probing harmonic oscillator, pulled 
out for clarity. 



motor can be used reversibly as a generator. As a gen- 
erator, the induced voltage is proportional to the angular 
velocity: e = a9. As a motor, the torque is proportional 
to the current supplied /: T = al. Notice that both re- 
lations involve the same proportionality coefficient, which 
depends on the physical parameters of the motor itself. 
The mechanical work produced by the motor against the 
granular gas per unit time is simply: w = T 6 = e I. 
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Fig. 2: The electrical sketch of the motor's command, and 
measurement set-up, that is to say the amplifiers-filters stage. 
The motor is represented by the induced voltage 'source' e and 
its internal resistance r. 

The motors command and measurement set-up are 
sketched in fig. [21 The harmonic oscillator is driven by 
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an AC current, supplied through a R = Ikil resistor by 
a sine voltage generator at fe = 13 Hz frequency. As R 
is large enough, the motor is driven by a periodic cur- 
rent (see time-series in fig. [3]). At such a low frequency, 
iron losses in the motor arc negligible. The voltages uq 
and ui are recorded by a 16 bits simultaneous acquisi- 
tion board at frequency fs = 1024 Hz. The amplifica- 
tors adjusts the level of the signals to that of the A/N 
converter, and anti-alias filters at 512 Hz. The results re- 
ported here arc obtained from one hom'-long recordings. 
The current I{t) = (uq — ui)/R and the induction voltage 
e(t) ~ ui — r I{t) are measured, r being the internal resis- 
tance of the motor (r ~ 5.8 fi). 

The power w{t) = T 9 ~ e I produced by the motor at 
13 Hz frequency is completely transferred to the granu- 
lafr gas. Although I{t) is a sine, e{t) fluctuates strongly 




Fig. 3: Top: time-series of the imposed sine current (/ ~ 
200/x^), and the induced voltage. For homogeneity and com- 
parison, both have been plotted in Volts: R I(t) (bold), and for 
readability, e has been enlarged by a factor 100 (thin). Bottom: 
in the same time-scale, the corresponding energy flux w{t) is 
plotted. The mean power given \s < w >~ 7.3 nW . The 
acceleration is a ~ 56 ms^^. 




f (Hz) 

Fig. 4: For the same as the previous measurement, the power 
spectral densities of e is plotted on the top, showing that the 
induced voltage (velocity) is broad-band noise. The excitation 
is visible at /e = 13 Hz. The bottom plot shows the power 
spectral density of the probing power w. 



ger multiples of the excitation period: r = n/ f^- Con- 
venient values are chosen as n = 10, 30, 50, 70, 90. Cor- 
responding histograms arc shown in fig. [S] They arc cal- 
culated over 50 bins. The wider histograms correspond 
to the smallest r. Starting from those histograms, eq. [5] 
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because of the beads collisions on the blade. The power 
w ^ el shows large fluctuations, negative whenever / and 
e have opposite sign (see fig. [3]). In such case, the reservoir 
is giving work to the blade-motor system. The fluctuations 
of power are widely distributed positively and negatively. 
Looking carefully at the time series, one cannot see in the 
fluctuations of the voltage any trace of the sine excita- 
tion. One can see in fig. H] the power spectral densities of 
the induced voltage e. It is broad band noise, reflecting 
that the momentum transfer from the granular gas is dis- 
ordered: no periodic contribution is visible except a small 
contribution at 13 Hz coming from the current. 

Results. — Histograms of the power Wr transmit- 
ted to the gas over several time-lags r are calculated: 
Writ) = - J e[t')I[t')dt' . Time-lags are chosen as inte- 



Fig. 5: Histograms of the power given to the gas for I ~ 200/iA, 
time-averaged over several values of r (D: 0.77 s, Q- 2.3 s, 



<1: 3.8s, o: 5.4s, [>: 6.9s.), i.e. 



to 8.9 Wtc, where 



Tc ~ 7.8 10" 
5.3 10"'-* W. 



w-r is normalised by the mean value < to >: 



is examined, plotting against uv the so-called asymmetry 
function: Mog f yi{-^ ) ) ■ Here and in the following, the 
natural logarithm only is considered. One can see in fig. |6] 
that these curves are approximately linear, at least for 
moderate Wr (t not too small), validating experimentally 
the prediction of FT. In some cases however, curves are 
observed to be slightly bent downward for large values of 
Wr. This is due to r being too small. 
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In order to test eq. [5J the slope of the asymmetry func- 
tion is evaluated. To avoid the fitting being dominated by 
extreme but less relevant values, it is performed for each 
T only over amplitudes Wr less than half of the maxima. 
A ratio slightly different, like 1/3 or 2/3, gives the same 
results. With 1/2, this procedure is stable for all exper- 
imental configurations considered here (change of /, and 
acceleration of the vessel) . 



FT takes the very simple form: 




Ec 



2 <Wr >' 



(3) 



WrI (w) 



Note that such a relation between the two sole moments 
of a Gaussian process is remarkable: the FT relates it by 
introducing the characteristic energy. 
In addition to the first method (slope of the asymmetry 
function), the calculation of Ec has been performed thanks 
to eq. [3] This second method to measure E^ valid under 
Gaussian hypothesis, does not require negative events. 
Those two methods are expected to give a close result, as 
Wr is close to Gaussian. Despite a large scattering of the 
slope measurements, fig. [7] shows that they are consistent 
within a range of about 20%, at least as far as negative 
events occur in large enough number to measure Ec- A 
discrepancy of this order of magnitude has been observed 
in all measurements. The process Wr is near Gaussian, 
but seems distinct at the present stage. 
Although the first method gives results more scattered 
than the second one, it will be preferred in the following 
as it is more general: no Gaussian assumption is needed. 



Fig. 6: The asymmetry function \log { yii™J, ) ) i® plotted 
against Wt/ < Wt >, for the same 5 time-lags as on previous 
figure (D: 0.77s, Q- 2.3s, <: 3.8s, o: 5.4s, [>: 6.9s, i.e. 980 
to 8.9 IQP Tc, where t^ ~ 7.8 10"'' s). (/ ~ 200^1^.) A linear 
fitting is performed, shown only after transient, for r ~ 2.3 s 
(deished line). 

The fitting parameter is the slope 1/ E^ where Ec is called 
characteristic energy. It is evaluated for different values of 
the time-lag r. It tends to converge toward a limit value 
at large r (see fig. Eland fig. [7]). 

For long times though, negative events become rare, caus- 
ing statistical uncertainty on Ec. For this reason, the 
parameters must be chosen such that < w > is small, 
and fiuctuation large. However, for short time-lags, as 
fiuctuations are dominant, this asymmetry-based method 
becomes unreliable. The whole procedure relies on a dcli- 

1 /9 

cate compromise. The fiuctuation rate < (Aw^) > / < 
Wr > is a decisive parameter for the reliability of this mea- 
surement. In this study, it is between 20 and 60. It is a 
very useful feature of this experiment to adjust < w,- > at 
will. This is labelled as the first method to measure Ec- 
It can be noticed on fig. [3] that the histograms are close to 
Gaussian, when r increases. This is expected in the limit 
of large r from central limit theorem for variables such as 
Wt, that result from a summation. It does not mean that 
w is Gaussian. Actually it is more like exponential, but 
the behavior at small t are not of interest here. Addition- 
ally, it has been checked that the variance goes as l/r as 
the mean is constant. 

For a Gaussian variable, it can easily be proved that the 
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Fig. 7: Ec measured at I ~ 2QQ^A, against t/tc, Tc being the 
correlation time of the process lii. The squares and circles refer 
to Ec measured according respectively to the first and second 
method, a — 56 ms^^. 

It is interesting to compute the correlation time Tc of w. 
The loss of correlation is rather fast, as the level of the 
periodic excitation is not much larger than the noise of 
the granular gas (see fig.lH top). In all the measurements 
shown, the correlation time Tc is constantly of the order 
of the ms, which is simply the inverse of the spectrum fre- 
quency span. This microscopic time Tc might simply be 
the mean fiight time of the beads between two collisions. 
The measurement described above is re-conducted for 
several values of the current /, at constant value of the 
vessel acceleration a = 41 ms"^. The curves shown in 
fig. [5] collapse to a unique value Ec — 1.7 lO^^J. Uncer- 
tainties are difficult to evaluate, but it is believed that the 
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dispersion of the curves is mainly due to the fitting process 
of the asymmetry function, and to statistical limitations. 
The collapse of the curves in fig. HI obtained at different 
/, can be interpreted as the following. Within a certain 
range of current amplitudes, Ec does not depend on the 
blade+motor system and its excitation. It is therefore a 
characteristic of the granular gas solely. 
In other words, this externally driven harmonic oscillator 
actually qualifies as a probe. 



2.5 



x10" 



"-^ 



1.5 



K] 



0.5 




2000 



4000 



6000 



8000 



Fig. 8: Ec against t/tc, for several values of the excitation 
current I (D : 1mA, <): 700 pA, and Q: 400 ^A), but the same 
acceleration: a = 41 ms^^. 

The range of 'acceptable' currents is however bounded. 
If / is too large, the fluctuation rate is small, angular 
velocity follows the torque, little negative events are ob- 
served: Ec can hardly be measured. If / is too small, the 
histograms are almost even, depreciating the resolution 
on Ec- For this reason, the time-lag on which Ec is 
constant and well defined is smaller for larger currents. 
This also can be seen in fig. [8] 

It has been shown so far that Ec is an energy 
characteristic of the granular gas, not of the probing, and 
how to measure it. No simple argument could be found 
to account for the actual numerical values of Ec- The 
following steps are to relate Ec to some parameters of the 
NESS granular gas itself. 



uncertainty, Ec = Ek- The characteristic energy of the 
granular gas may be interpreted as the kinetic energy of 
the probe coupled to the granular gas. 
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Fig. 9: Ec is plotted against the kinetic energy of the 
blade-frotor system: Eh = ^M < (P >. Dashed line is the 
best linear fit through 0. The acceleration of the vessel is 
varied between 41 ms^^ and 59ms~^, at constant frequency 
/^xc. = 40 Hz, I = 200 ^A unchanged. 

Another parameter that can easily be varied is the ac- 
celeration a imposed by the shaker to the beads. How can, 
dimensionally, an energy be related to an acceleration ? 
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Varying the power supply of the shaker Pe^c., one 
easily vary the acceleration of the vessel and therefore the 
agitation of the granular gas. A calibration of the motor 
is performed, giving the coefficient a = e/9 ~ 0.0029 
with a good accuracy. The calibration of the moment of 
inertia gives M = 3.2 10~^kgm^, with an uncertainty 
of about 15%. The kinetic energy of the blade-|-rotor 
system E^ = ^A/ < 9^ > can be calculated, from the 
variance of the induced voltage < e^ >, recorded with no 
excitation. In fig. [SJ Ec is plotted against Ek for several 
values of the acceleration. A linear fitting, forced at 0, 
gives a proportionality coefficient of 0.89. Within the 



Fig. 10: Ec against (a/(27r/exc.))^, for several values of the 
acceleration (between 41 ms^^ and 59 ms^^). The frequency 
/^xc. = 40 Hz, / = 200 ^A unchanged. 

Supposing some sort of thermalisation between the probe 
and the granular gas, Ec can be written tentatively in 
the form: Ec = E - Eq, with E =< mv'^ >. The 
mass m accounts for a certain mass of beads, and v their 
velocity. The brackets represent an average on an un- 
specified number of particles in the vicinity of the blade. 
Eq is the minimal energy needed to take a large enough 
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number of beads up to the probe. For E < Eq, there 
is no Ec measurable, simply because no beads hit the 
blade. Still dimensionally, this relation can be rewritten 
as: Ec = ?7i'(a/(27r/cxc.))^ — Eq. No negative Ec nor Eq is 
to be considered. Apart from small values, data arc con- 
sistent with this naive picture, at least at the present level 
of uncertainty, as can be seen in fig. 1101 Uncertainties on 
Ec arc hard to quantify, but may be large. In these con- 
ditions it is vain to further investigate Eq and m' , which 
are system-dependent. 

Conclusion and perspective. The system pre- 

sented here is very simple. Nevertheless, it allows very 
sensitive measurements and gives access to important and 
topical issues in non-equilibrium statistical mechanics. 
It consists of a driven harmonic oscillator, coupled as a 
probe to a granular gas. It is composed by a DC-motor 
with a blade and a torsion spring on its axis. The gran- 
ular gas is maintained in a gaseous state by an external 
power supply Pcxc. from a shaker. It has an extremely 
large number of degrees of freedom, which guaranties an 
efficient chaotic loss of memory. It forms a NESS bath, 
to which the harmonic oscillator thermalises. The gran- 
ular gas mimics quite well a real gas, yet with important 
differences: the number of particles is not extremely large 
(fluctuations are important), the collisions are dissipative 
(a power supply is necessary). 

The DC-motor driving the harmonic oscillator is used as 
actuator and sensor altogether, to probe the granular gas. 
The electro-mechanical relations in a DC-motor arc such 
that no calibration is needed to measure the power deliv- 
ered to the bath. This property is fortunate, as calibration 
is a source of errors. 

Through this device, a small energy flux w is injected into 
the granular gas at controlled torque. Note that < w > is 
of the order of a few 10~^ W, i.e. resolution is much less. 
It is extremely small in such a macroscopic experiment ! 
The fluctuations of this perturbation power is shown to 
likely verify the Stationary State Fluctuation Theorem, 
although the heat bath is non-thermal. However, it must 
be noted that this study does not strictly rely on this the- 
orem: the asymmetry function is used to define the output 
quantity Ec- The characteristic energy does not depend 
on the imposed torque amplitude within a certain range. 
One concludes that Ec is characteristic of the granular gas 
itself, i.e. this small system is actually a probe. 
For different values of the power Pexc. supplying the NESS, 
the characteristic energy Ec is measured, and found equal 
to the kinetic energy of the swivel device. The dependance 
in the acceleration is examined. Despite important scat- 
tering of the data, a clear tendency to increase is visible, 
compatible with a quadratic dependance in a. This cor- 
roborates the intuition. 

A quantitative characterisation of this relation definitely 
requires improvements of the experimental set-up and the 
measurement procedure to measure the asymptotic Ec. 
Considering the FT on a small perturbative energy flux 



instead of that supplying the NESS is a biased opinion 
adopted here, which appears fruitful. It allows in prin- 
ciple to measure negative fluctuations of injected power 
with arbitrarily long statistics, for a granular gas as far as 
desired from equilibrium. This key point was discussed by 
Zamponi [7] as a limiting benchmark for experiments. 
One interesting question is whether Ec obtained via the 
FT in the present study is the same or different as that 
obtained using the fluctuations of the power supply Pcxc, 
that maintains the gas in its NESS. (Like in most of pre- 
vious studies, for instance [5]). 

It would be interesting to clarify the physical significance 
of Ec. The characteristic energy is often interpreted, in the 
context of FT, in terms of an effective (non-equilibrium) 
temperature: Ec = fceToff.. Ec has been related to the 
kinetic energy of the probe, but experimental evidence of 
a relation with the kinetic energy of the gas in the vicinity 
is still needed to justify thermalisation. Such a measure- 
ment is definitely needed to resolve the issue. But such 
study would require different technics, like PIV. 

Overextending the analogy with a gas, it is worth notic- 
ing that the relation obtained in fig. IH] differs from the pre- 
diction of the equipartition theorem, by a factor 2. 
In the present experiment, such interpretations deserve 
caution, as the forcing is non-Gaussian (for short times), 
the granular gas is intrinsically dissipative and has an irre- 
versible dynamics (the heat bath is non-thermal). Rigorous 
interpretation of these observations remains to be done. 
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